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Abstract 

We study the phase transition of thin films in the three-dimensional XY 
universality class. To this end, we perform a Monte Carlo study of the 
improved two-component 0^ model, the improved dynamically diluted XY 
model and the standard XY model on the simple cubic lattice. We study 
films of a thickness up to Lq = 32 lattice spacings. In the short direction of 
the lattice free boundary conditions are employed. Using a finite size scaling 
(FSS) method, proposed recently, we determine the transition temperature 
with high accuracy. The effectively two-dimensional finite size scaling be- 
haviour of the Binder cumulant C/4, the second moment correlation length 
over the lattice size 5,2nd/ L, the ratio of the partition functions with anti- 
periodic and periodic boundary conditions Za/Zp and the helicity modulus 
T clearly confirm the Kosterlitz-Thouless nature of the transition. We anal- 
yse the scaling of the transition temperature with the thickness Lq of the 
film. The predictions of the renormalization group (RG) theory are con- 
firmed. We compute the universal ratio of the thickness of the film Lq and 
the transversal correlation length in the three-dimensional thermodynamic 
limit at the Kosterlitz-Thouless transition temperature of a film of thickness 
Lq: [Lq^kt/S,t\* = 1.595(7). This results can be compared with experimental 
results on thin films of *^He near the A-transition. 
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1 Introduction 



In the neighbourhood of a second order phase transition the correlation length 
diverges as 

e^/±ir% (1) 

where t = {T — Tc)/Tc is the reduced temperature, /+ and /_ are the amplitudes in 
the high and the low temperature phase, respectively, and u is the critical exponent 
of the correlation length. In the neighbourhood of the transition also the behaviour 
of other quantities is governed by power laws. E.g. the specific heat behaves as 

(:7~A±|t|-" + S, (2) 

where B is an analytic background, which has to be taken into account here, since 
the critical exponent a of the specific heat is negative for the three-dimensional XY 
universality class, as we shall see below. Critical exponents like u and a and ratios 
of amphtudes such as f+/f- and A^/A_ assume universal values. I.e. they are 
supposed to assume exactly the same value for any system in a given universality 
class. This can be understood in the framework of the Renormalization Group 
(RG). A universality class is characterised by the dimension of the system, the 
range of the interaction and the symmetry of the order parameter. For reviews on 
critical phenomena and the Renormalization Group see e.g. [U [2l [3l H] . 

At temperatures below the A-transition, ^He becomes superfluid. The order 
parameter of the A-transition is the phase of a wave function. Therefore it is sup- 
posed to share the three-dimensional XY universality class, which is characterized 
by the 0(2), or equivalently U(l), symmetry of the order parameter. The most 
accurate experimental results for critical exponents and amplitude ratios are ob- 
tained for the A-transition of ^He. E.g. under the condition of micro-gravity, the 
specific heat has been studied for reduced temperatures as small as t ~ 5 x 10~^° 
[SI E], resulting in a = —0.0127(3), corresponding to u = 0.6709(1). Note that 
the exponents of the specific heat and the correlation length are related via the 
hyperscaling relation du = 2 — a, where d is the dimension of the system. Ex- 
periments on earth, measuring the specific heat and the second sound to obtain 
the superfluid density have resulted in accurate estimates of the exponent of the 
correlation length u = 0.6717(4) and u = 0.6705(6) in refs. OE], respectively. For 
a recent review and an outlook on future experiments in space-stations, see ref. [9]. 
The experimental results for critical exponents are in reasonable agreement with 
the most precise theoretical prediction for the three-dimensional XY universality 
class: u = 0.6717(1) obtained from a combined Monte Carlo and high temperature 
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series study of improved lattice models [TO]. For a summary of theoretical results 
obtained with various methods (e.g. field-theoretic methods) see ref. [3]. Also in 
the case of universal amplitude ratios, there is reasonable match between theoretical 
results for the three-dimensional XY universality class and experimental studies of 
the A-transition. See refs. [Hj [12] and refs. therein. 

The discussion above refers to the thermodynamic limit. It is a natural question 
how the behaviour at the transition is altered by a finite extension of the system. It 
has been addressed both experimentally and theoretically. In systems with a finite 
extension in all directions, there can not be any singularity such as eqs. ( TT1I21) . As 
a remnant of these singularities there remains a peak in the neighbourhood of the 
transition. With increasing linear extension the hight of the peak increases and the 
temperature of the maximum approaches the critical temperature. 

The situation might be different, if the extension in some of the directions stays 
infinite. In the case of one infinite direction and two finite ones the system becomes 
effectively one-dimensional. For the type of interactions that we are dealing with 
here, there is no phase transition at a finite temperature for a one-dimensional 
system. Hence one expects that the behaviour of the effectively one-dimensional 
system is qualitatively the same as that of a system which is finite in all directions. 

In the case of a thin film geometry, i.e. one direction with finite extension and 
two infinite ones, we expect that the system becomes effectively two-dimensional in 
the neighbourhood of the transition. Hence, there will be still a phase transition, 
however it belongs to the two-dimensional universality class. I.e. in our case of U(l) 
symmetry of the order parameter, a Kosterlitz-Thouless (KT) transition [131 El IIS] 
is expected. 

The behaviour of finite systems in the neighbourhood of a continuous transition 
is described by finite size scaling (FSS). For a reviews see [TBIIIZ]. In essence, close 
to the transition the bulk correlation length ^, i.e. the correlation length in the 
three-dimensional thermodynamic limit, and the extension of the finite system are 
the only relevant length scales. I.e. the physics of finite systems is governed by 
the ratio Lq/^. In particular, one expects that, independent on the thickness Lq of 
the thin film, the effectively two-dimensional phase transition occurs at a universal 
value of Lq/^. One should note that this universal value depends on the type of 
boundary conditions that are applied in the finite direction of the system. Using 
eq. ([ID it follows that [TH HQ] 



where Tc^s^ is the transition temperature of the three-dimensional system and Tkt{Lq) 
the temperature of the KT transition of the thin film of thickness Lq. 




(3) 
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The predictions of finite size scaling have been checked by various exphcit cal- 
culations and experiments fl6[ [T7] . For a recent review of experimental results near 
the A-transition of ^He and ^He-^He mixtures see [20]. For thin films in the three- 
dimensional XY universality class a multitude of field theoretic calculations have 
been performed which allow for comparison with experiment [21]. Some important 
aspects of the behaviour of thin films are however not accessible using field theoretic 
methods. Among these is the KT transition, which is the focus of the present work. 

In thin films of '^He on various substrates, the order parameter is vanishing at 
the boundary of the film. This has to be taken into account in the Monte Carlo 
simulation of lattice models with film geometry. The simplest way to obtain a 
vanishing field at the boundary are so called free boundary conditions. For a precise 
definition see the section below. An alternative are so called staggered boundary 
conditions. The numerical results of ref. p2] confirm that free and staggered 
boundary conditions lead the same results for universal quantities. 

On top of the restricted geometry, free boundary conditions might introduce 
new physical effects; For reviews on surface critical phenomena see refs. [231 |2l]. In 
fact, free boundary conditions lead to additional corrections to scaling. The leading 
one is oc Lq^ [2^; it can be cast in the form Lo^ff = Lq + Lg. These corrections 
come in addition to those oc Lq'^ , oc Lq'^ , which are predicted for finite systems 
in general, irrespective of the type of the boundary conditions [16]. Note that the 
numerical value of the correction exponent is = 0.785(20) [TD]; similar results are 
obtained with field-theoretic methods; see e.g. the review [1]. The information on 
uj' is rather sparse; following ref. [26] uj' 2uj. Fitting Monte Carlo data, it might 
be difficult to disentangle corrections oc Lq*^ and oc Lq ^, leading to sizable errors in 
the extrapolation Lq — oo. In order to avoid this problem, and to clearly show the 
existence of corrections oc Lq ^ due to the free boundary conditions, we have studied 
improved models. In these models the amplitude of corrections oc Lq vanishes, or 
in practise, it is so small that its effect can be ignored. The precise definition of 
these models is given below. 

In the literature on can find only a few Monte Carlo studies of lattice models 
that address the scaling of the KT temperature with the thickness of the film. Janke 
and Nather [27] have studied the Villain model on the simple cubic lattice with free 
boundary conditions in the short direction. They have determined the temperature 
of the KT transition from the behaviour of the correlation length and the magnetic 
susceptibility in the high temperature phase of the thin films. They have studied 
films of a thickness up to Lq = 16. Fitting the transition temperature with eq. ([2D 
they find a value for the exponent v that is about 5% too large compared with the 
estimates of v discussed above. 
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Schultka and Manousakis [28] have studied the standard XY model with periodic 
boundary conditions in all directions. They studied lattices up to a thickness Lq = 
10. They determined the KT transition temperature using the effectively two- 
dimensional finite size scaling behaviour of the helicity modulus. They conclude 
that their results are consistent with eq. (I3l), using z/ = 0.6705, which was the best 
experimental estimate of the correlation length exponent at the time. The same 
authors [291 EH] have studied the standard XY model with staggered boundary 
conditions in the short direction. Also in this case, they have determined the 
KT transition temperature using the effectively two-dimensional finite size scaling 
behaviour of the helicity modulus. They studied films of the thicknesses Lq = 
4, 8, 12, 16 and 20. They find that their data can be fitted with eq. ([3]) using u = 
0.6705. However in this case, that requires that Lq in eq. ([3]) is replaced by Lo^^ff = 
Lq + Lg, with Lg = 5.79(50). They have also reanalysed the data of ref. [27] this 
way. They find that also these results are compatible with u = 0.6705, once Lq is 
replaced by -Z^o,e// with Lg = 1.05(2). One should note however that Schultka and 
Manousakis do not take into account corrections oc Lq^ which should be present in 
the standard XY as well as in the Villain model. 

There are further Monte Carlo studies of thin films using lattice models in 
the three-dimensional XY universality class. These studies focus on the specific 
heat [22] and refs. therein, the thermal resistivity [31] and refs. therein, or the 
thermodynamic Casimir force [321 [33]. Throughout, the thickness of the films that 
had been studied is less or equal to Lq = 24. 

This paper is organized as follows: In the following section, we define the lattice 
models that we have studied. Then we introduce the observables that we have 
measured. In section \5\ we study the effect of free boundary conditions in finite size 
scaling (FSS) directly at the critical temperature of the three dimensional system. 
In section [6] we briefly summarize the method proposed in ref. \3^ to determine 
the KT transition temperature. Then we determine the KT transition temperature 
for a large range of the thickness of the film. We demonstrate that the transition 
of the thin films is indeed of KT nature. To this end we compare the behaviour 
of various phenomenological couplings, also called dimensionless quantities in the 
following, at the transition temperature with predictions from KT theory and from 
numerical results for two-dimensional systems that are known to undergo a KT 
transition. Next we study the scaling behaviour of the KT transition temperature 
with the thickness of the film. Then we compute universal amplitude ratios that 
relate the thickness of the film with the correlation length of the bulk system at 
the temperature of the KT transition of the thin film. We compare our estimate 
for the universal amplitude ratio with that of experiments on thin films of ^He. 
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2 The models 



We study various models with 0(2)-symmetry on a simple cubic lattice. We consider 
systems with film geometry. We shall label the sites of the lattice by x = {xq, xi,X2)- 
The components of x might assume the values Xj G {1, 2, Lj}. We simulate 
lattices of the size Li = L2 = L and Lq <^ L. In 1 and 2-direction we employ 
periodic boundary conditions and free boundary conditions in 0-direction. This 
means that the sites with xq = 1 and xq = Lq have only five nearest neighbours. 
This type of boundary condition could also be interpreted as Dirichlet boundary 
conditions with as value of the field at Xq = and Xq = i^o + 1- 
The standard XY model is given by the Hamiltonian 

HxY = -13 ^ s^Sy , (4) 

<x,y> 

where is a unit-vector in M?. < x,y > denotes a pair of nearest neighbour sites 
on the lattice. In our convention, the inverse temperature (3 is absorbed into the 
Hamiltonian. The Boltzmann factor is given by exp(— Tixy)- The best estimates 
of the inverse transition temperature that are quoted in the literature are f5c = 
0.454165(4), 0.454167(4), 0.4541659(10) and 0.4541652(11) in refs. [35l [Ml EH [10] , 
respectively. In the following we shall assume (5c = 0.4541655(10) which is, roughly, 
the average of the estimates given by refs. [371 [TO] - 

A generalization of the XY model is the 0^ model on the lattice. The Hamilto- 
nian is given by 

n^. = -f3j2$.-$y + J2[^' + - 1)'] ' (5) 

<x,y> X 

where the field variable (px is a vector with two real components. The partition 
function is given by 



n 



exp(-7^04). (6) 



In the limit A — >■ 00 the field is fixed to unit length; i.e. the XY model is recovered. 
For A = we get the exactly solvable Gaussian model. For < A < 00 the 
model undergoes a second order transition that belongs to the XY universality 
class. For a discussion see e.g. [S^. Numerically, using Monte Carlo simulations 
and high-temperature series expansions, it has been shown that there is a value 
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A* > 0, where leading corrections to scaling vanish. Numerical estimates of A* 
given in the literature are A* = 2.10(6) [39], A* = 2.07(5) [IQ] and most recently 
A* = 2.15(5) [To]. The inverse of the transition temperature has been determined 
accurately for several values of A using finite size scaling (FSS) |LlQj. We shall 
perform our simulations at A = 2.1, since for this value of A comprehensive Monte 
Carlo studies of the three dimensional system in the low and the high temperature 
phase have been performed [T0 | ITT j [12]. The inverse critical temperature at A = 2.1 
is Pc = 0.5091503(6) [10]. Since A = 2.1 is not exactly equal to A*, there are still 
corrections oc L~'^, although with a small amplitude. In fact, following ref. [ID], it 
should be by at least a factor 20 smaller than for the standard XY model. 
The dynamically diluted XY model [iQl [10] is given by 

HddXY = -/5 ^ 0. ■ 0, - D ^ (7) 

(xy) X 

with the local measure 

rf/x(0,.) = #(2) 

i.e. \(j)x\ is either or 1, and the partition function 

j nc?/^(0x) exp(-7YddXY) • (9) 



^ddXY 



In the limit oo the standard XY model is recovered. For D < Dtri the model 
undergoes a first order transition and for D > Dfri a second order transition in 
the XY universality class. A mean-field calculation gives Dtri = 0, while from a 
improved mean- field calculation we get Dtri < [ID]- Therefore it is quite save to 
assume that indeed Dtri < for the ddXY model. 

From Monte Carlo simulations and high temperature series expansion we know 
that there is a Dtri < D* < (yo such that leading corrections to scaling vanish. 
Numerical estimates are D* = 1.02(3) [40j and D* = 1.06(2) [lU]. Here we shall 
present some preliminary results for D = 1.02. At this value of D, the inverse critical 
temperature is Pc = 0.5637963(4). Also here, the amplitude of leading corrections 
to scaling should be by at least a factor 20 smaller than for the standard XY model 

m. 
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3 The observables 



The total magnetisation is defined by 



(10) 



Note that in the case of the XY model we have to replace (f)x by Sx here and in the 
following definitions. The magnetic susceptibility in the high temperature phase, 
for vanishing external field, is given by 



X 



L0L1L2 



where (...) denotes the expectation value with respect to the Boltzmann factors 
defined in the previous section. The Binder cumulant is defined by 



We also consider the generahzation 



{(Mr) 

{MY 



(M2)3 



(12) 



(13) 



We have computed the second moment correlation length for the 1- and 2- 
direction (i.e. the large ones). The second moment correlation length in 1-direction 
is defined by 



s2nd,l 



4 sin(7r/Li)^ 



where 



1 

L0L1L2 



E 



exp I t 



271X1 



0a 



(14) 



(15) 



is the Fourier transform of the correlation function at the lowest non-zero momen- 
tum in 1-direction. In our simulations, we also have measured F2 in order to reduce 
the statistical error. 

The helicity modulus T gives the reaction of the system under a torsion. To 
define the helicity modulus we consider a system, where rotated boundary conditions 
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are introduced in one direction: For xi = Li and yi = 1 the term (j)x4>y in the 
Hamiltonian is replaced by 

0x ■ Ra(py = (cos(a)0W + sm{a)(f)^^^) + 01') (- sin(a)0« + cos(a)0(2)) . (16) 

The hehcity modulus is then given by 

Li d"^ log Z{a) 



T 



L0L2 da^ 



(17) 



Q = 




Note that we have skipped a factor of T compared with the standard definition 
Defined this way, the helicity modulus has the dimension of an inverse length. In 
the literature = 1/T is referred to as transversal correlation length. The helicity 
modulus can be directly evaluated in the Monte Carlo simulation. Following eq. (3) 
of ref. [H] one gets for the models considered here 

T = ^ I 44+(0.1.0) /- r f r 

1^0-^1-^2 \ / -t^O-^1-^2 

(18) 

In addition we have measured the analogous quantity for the 2-direction. A problem 
of the estimator f|T8l) of T is that for a fixed number of measurements and LjT 
fixed, the statistical error is increasing as the critical point is approached; see e.g. 
ref. [12]. In section H] below we shall discuss how this problem can be reduced 
to some extend. From the definition f[T7l) one reads off that LqT is dimensionless. 
Therefore, in the case of thin films the behaviour of LqT has to be compared with 
the prediction for two-dimensional systems in the neighbourhood of the Kosterlitz- 
Thouless transition. 

The ratio of partition functions Za/Zp is a useful quantity in the numerical study 
of phase transitions. Here Zp is the partition function of a system with periodic 
boundary conditions in two directions of the lattice and free boundary conditions 
in the remaining direction, while Za is the partition function of a system with 
anti-periodic boundary conditions in one of the directions and periodic and free 
boundary conditions in the other two directions. Anti-periodic means that, e.g. 
in the 1-direction, at the boundary, i.e. for xi = Li and yi = 1 the term (I)x4'y 
is replaced by —4>x4>y in the Hamiltonian. This can also be viewed as a rotation 
by TT. The ratio of partition functions can be efficiently measured in Monte Carlo 
simulations using a special version of the cluster algorithm |13] . The implementation 
of the algorithm that is used here follows the discussion given in Appendix A 2. of 
ref. sol. 
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4 The Monte Carlo Algorithm 



For our simulations we use a composition of local and cluster algorithms. The 
implementation is similar to that used in ref. [ID]. See Appendix A 1. of ref. 
[IQ] for details. Applied to the standard XY model cluster algorithms are ergodic. 
However, in the case of the ddXY model and the (p'^ model additional local updates 
are needed to change \(f)x\- As cluster algorithms, we have used both the single- 
cluster IH] and the wall-cluster algorithm [45j . 

Let us briefly discuss the details of the local updates for the case of the cp'^ model 
for which most of our simulations were performed. In the case of the Metropolis 
update a proposal is generated as 

05 = (19) 

with s = 3 (the choice of s is such that very roughly the acceptance rate is 50%) 
and 

C = 5^ 4 > (20) 

y.nn.x 

where y.nn.x means that y is a nearest neighbour of x. The idea of this proposal 
is to change also the second component without using a random number. The 
probability to accept this proposal is Pace = min[l, exp(— Aif)]. This update is 
followed by a second one at the same site, where the role of the two components is 
exchanged. Using these updates we run through the lattice in typewriter fashion. 
Going through the lattice once is referred to as one sweep in the following. 

We have also implemented overrelaxation updates for all models that we have 
simulated. These are given by 

0; = 0. - 2^^$. , (21) 

X 

where 

y.nn.x 

Note that these updates do not change the value of the Hamiltonian and therefore 
no accept/reject step is needed. Hence it is computationally quite cheap since no 
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random number and no evaluation of exp() is needed. Also in the case of the over- 
relaxation update we run through the lattice in typewriter fashion. As the cluster 
updates, the overrelaxation update does not change \(t)x\- The main motivation to 
use the overrelaxation update is that it allows to reduce the variance problem of the 
helicity modulus to some extend. Since the overrelaxation update is quite cheap in 
terms of CPU-time, we performed several sweeps with the overrelaxation update in 
each update cycle. After each of these sweeps the second term on the right hand 
side of eq. f fTSj) is measured. For a lack of human time we did not carefully tune 
the number of overrelaxation sweeps per update cycle. 

As random number generator we have used the SIMD-oriented Fast Mersenne 
Twister (SFMT) |16] generator. In particular, we use the genrandjres53() function 
that produces double-precision output. 

5 Finite size scaling at pc,3d 

First we performed simulations at the transition temperature of the three-dimensional 
system. This way we avoid possible difficulties related with the KT transition and 
hence can see more clearly the effects caused by free boundary conditions. Correc- 
tions caused by the free boundaries should qualitatively not depend on Li and L2, as 
long as Li and L2 scale with Lq. Therefore we consider lattices with Lq ~ = L2, 
allowing us to study a relatively large range in Lq. 

We have simulated the 0^ model at A = 2.1 at the best estimate of the inverse of 
the critical temperature of the three-dimensional system: f3 = 0.5091503 [lOj. We 
consider the lattice sizes Li = L2 = 8, 12, 16, 24, 32 and 48. In addition to Lq = Li 
we have simulated Lq = Li — 1. Throughout we have performed 10^ measurements. 
For each of these measurements, one boundary flip update, several single cluster 
updates, one Metropolis sweep, and three overrelaxation sweeps were performed. 
Note that this composition is essentially chosen ad hoc. In total, these simulations 
took about 6 month of CPU time on one core of a Xeon(R) CPU 5160 running at 
3 GHz. We have measured the dimensionless quantities Za/Zp, ^2nd/L, U4, and Uq 
as defined in section [31 For all these quantities we have determined the coefficients 
of the Taylor series in /3 around /5 = 0.5091503 up to third order. This allows us to 
evaluate these quantities in the neighbourhood of (3 = 0.5091503. Our results for 
Za/Zp and ^2nd/L are plotted in figure [H The results for U4 and Uq are given in 
figure O In the case of Lq = Li corrections are clearly visible for all four quantities 
considered. Plotted as a function of l/Li, the data for Lq = Li follow roughly a 
straight line, indicating that the corrections are proportional to L^^, as expected. 
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Figure 1: The 0^ model at A = 2.1 ttid (3 = 0.5091503. The upper (lower) 
figure gives Za/Zp {^2nd/L) as a function of 1/Li. We have simulated systems with 
Lq = Li = L2 and Lq + 1 = Li = L2. 
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Figure 2: Same as previous figure, but for f/4 and Uq. 



Table 1: Fits of Za/Zp with the ansatz ( f23j) . Li^rnin is the minimal lattice size 
included into the Gt. For a discussion see the text. 





{zjZpT 


Ls 


xVd.o.f. 


8 


0.70474(4) 


1.117(1) 


111.45 


12 


0.70364(5) 


1.068(2) 


11.84 


16 


0.70326(7) 


1.046(3) 


1.68 


24 


0.70311(12) 


1.033(7) 


0.18 



Table 2: Fits of C,2nd/L with the ansatz / f23|) . Li^min is the minimal lattice size 
included into the fit. For a discussion see the text. 



Ll,min 




Ls 


xVd.o.f. 


8 


0.36988(2) 


1.066(1) 


16.97 


12 


0.36988(3) 


1.058(2) 


4.98 


16 


0.36988(4) 


1.053(4) 


2.80 


24 


0.36986(7) 


1.049(8) 


0.51 



In addition, we give results for lattices with Lq = Li — \. Here, for the largest 
values of Li, the curves are almost flat for all quantities considered. E.g. the plots 
suggest that the dominant corrections can be explained by an effective thickness 
Lo,e// = Lq + Ls with Ls ^ 1. 

In the following we put this discussion on a more quantitative level by performing 
a sequence of fits of our data for the dimensionless quantities. 

First we have fitted the dimensionless quantities with the ansatz 

R{Li,Li-L^)=R* + c{Li-Lo)L^^ , (23) 

where R* , c(0) and c(l) are the parameters of the fit. Lg is given by the zero of 
c(Li — Lo). Using our data, we have linearly interpolated/extrapolated c(0) and 
c(l) resulting in 

L. = . (24) 

c(0)-c(l) 

Results of such fits are given in the tables [H [21 [3] and |4] for Za/Zp, i^nd/ 
and Uq, respectively. 
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Table 3: Fits ofU^ with the ansatz ( f23j) . Li^rnin is the minimal lattice size included 
into the fit. For a discussion see the text. 



Ll,min 


u: 


Ls 


xVd.o.f. 


8 


1.53578(7) 


1.188(5) 


28.39 


12 


1.53545(9) 


1.119(7) 


7.57 


16 


1.53540(12) 


1.097(12) 


6.49 


24 


1.53531(20) 


1.055(24) 


0.31 



Table 4: Fits of Uq with the ansatz l[23\) . Li,mm is the minimal lattice size included 
into the fit. For a discussion see the text. 



Ll,rain 


TJ* 


Ls 


xVd.o.f. 


8 


2.9573(3) 


1.224(5) 


33.83 


12 


2.9556(5) 


1.145(8) 


8.57 


16 


2.9553(6) 


1.119(12) 


6.84 


24 


2.9547(10) 


1.070(24) 


0.63 
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In all four cases, an acceptable x^/d-o.f. is only reached for Li mm = 24. From 
the fits with Li^rnin = 24 we get quite consistent results for for the four quantities 
considered. They range from Lg = 1.033(7) for Za/Zp up to Ls = 1.070(24) for Uq. 

In order to check the possible error due to the uncertainty of (3c, we have repeated 
the fits for Li,™„ = 24 for (3 = 0.5091509. The result for e.g. Za/Zp is (Za/Zp)* = 
0.70294(12) and Lg = 1.0277(7). I.e. the effect is of similar size as the statistical 
error. 

Next we performed fits with the ansatz 

i?(Li, Li - Lo) =R* + Ci(Li - Lo)Li' + c^iL, - 1^)1^ , (25) 

where now R*, Ci(0), Ci(l), C2(0) and C2(l) are the free parameters of the fit. Here 
we have not included explicitly L~'^ corrections, since at the level of our numerical 
resolution u' is rather close to 2. 

Using this ansatz, we get x^/d-o.f. ~ 1 already for Li^rnin = 8 for all four 
quantities that we study. The results for Lg from these fits are slightly smaller 
than those from the ansatz ([25]): = 0.955(5), 1.040(6), 0.991(19) and 0.998(18) 
from Za/Zp, ^2nd/L, f/4 and Uq, respectively. In the case of Za/Zp, where the result 
for Lg is the smallest, we get Lg = 0.969(11) from Li^rnin = 12. I.e. it is moving 
towards the other results. Also here, we have checked the effect of the uncertainty 
in pc. E.g. for (3^ = 0.5091509 and Li,„i„ = 8 one obtains Lg = 0.950(5), 1.035(7), 
0.981(18), and 0.988(18), from Za/Zp, ^2nd/L, U4 and Uq, respectively. I.e. the 
effect is of similar size as the statistical error at fixed f3. 

As our final result we quote 

Lg = 1.02(7) (26) 

where the error-bar covers all the results of the fits reported above. 

Here we did not explicitly check the possible error due to residual corrections 
oc at A = 2.1. However, corrections oc L~^ should affect different phenomeno- 
logical couplings in a different way. Therefore the variation of Lg obtained from dif- 
ferent quantities should also provide an estimate of the error resulting from residual 
corrections oc L^'^. 

Our final values for the fixed point values of the dimensionless quantities are 
(Za/Zp)* = 0.7030(5), {^2nd/L)* = 0.36986(10), = 1.5352(2), and = 2.954(1). 
These values can be compared with those of a L^ system with periodic bound- 
ary conditions in all three directions [lOj: (Za/Zp)* = 0.3203(1)[3], (^2nd/L)* = 
0.5924(1)[3], Ul = 1.2431(1)[1], and = 1.7509(2)[7], where the number in [] gives 
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the systematic error due to scaling corrections. The effect of the boundary condi- 
tions is very well visible. In the case of {Za/Zp)* the value even differs by more 
than a factor of two. 

The numbers quoted here, could be used to check whether other boundary con- 
ditions, e.g. staggered boundary conditions, are equivalent with free boundary 
conditions. One would expect a different value of Lg but the same values of R* for 
equivalent boundary conditions. 



6 Finite size scaling at the KT transition 

Finite size scaling is an efficient method to locate and to verify the nature of the KT 
transition in two-dimensional systems. To this end, usually the helicity modulus T 
is studied. Following ref. jlTj it behaves as 

2 1 1 

TT TT (InL + C) 

at the transition temperature. Recently, we have pointed out [48j that for lattices 
with periodic boundary conditions there are numerically small contributions from 
winding configurations that are not taken into account in eq. (I27j) . Including these 
contributions one gets [IE] 



_ 0.318899454... , , 

^LMransition = 0.63650817819... + ^ ^ ^ + ... (28) 

for an lattice with periodic boundary conditions in both directions at the KT 
transition. In the appendix we derive for the ratio of partition functions 

iZJZp)L2,trans^uon = 0.0864272337... - ^'^^^]^^^"' + - (29) 

for an lattice at the KT transition. We have also determined the leading finite 
size behaviour of the second moment correlation length over the lattice size 



2nd 



0.212430... , , 

= 0.7506912... + —— — Tr+ - (30) 

L"^ transition ^ ~^ 



and more recently for the Binder cumulant 



, , 0.017922(5) , , 

Ui,L\transition = 1.018192(6) - ^ ^ ^ + - ■ (31) 
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Eqs. ( !28| [30!) describe quite well the behaviour of numerical data for T and ^2nd/ L 
for various models down to rather small lattice sizes. In contrast to that, on the 
accessible lattice sizes, f/4 is decreasing with increasing L, while eq. (I3TI) suggests 
that it should increase. This behaviour is explained by disorder caused by vortex 
pairs at a distance of the order L/2. Therefore a term C2/(lnL + C)^ has been 
added to eq. fl3T|) . The coefficient C2 and C have been determined from fits to Monte 
Carlo data. The result for C2 is consistent for different models. This confirms the 
assumption that also the subleading correction in the Binder cumulant is universal 



In ref. [3l] we suggest that the KT transition of a model can be accurately 
determined by matching ^2nd/ L and the Binder cumulant with that obtained 
from a model where (3kt is accurately known. In ref. [31] we have studied most 
accurately the dual of the absolute value solid-on-solid (ASOS) model. In this 
model, the KT transition occurs at /3 = 0.80608(2) [49j. 

The data of the dual of the ASOS model for L > 96 are well described by 

(3 = 0.80608 : 

. ..o..^ 0.017922 0.06769 
U^ASosiL) = 1.018192 -— — - + 



lnL-1.18 (InL- 1.18)2 

^2nd,Asos{L) ^ 0.7506912 + ^:^1^^. (32) 
L In L + 0.573 ^ ' 

To check the error due to the uncertainty of Pkt one should repeat the matching 
with 

/3 = 0.80606 : 

.n.o.no 0.017922 0.06730 
U4,AsosiL) = 1.018192 -— — — + 



InL- 1.193 (InL- 1.193)2 

^2ndAsos{L) 0.212430 , , 

^2nd,AbusK J ^ 0.7506912 + — . 33 

L lnL + 0.557 ^ ' 

In order to match the thin films with these expressions one has to adjust (3 and a 
scale factor for the lattice size. 



7 KT transition in thin films 

Here we have studied the KT transition of thin films of the (j)^ model at A = 2.1 using 
the matching method discussed above. We have simulated lattices of the thickness 
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Lq = 4, 6, 8, 12, 16, 24 and 32. In order to get the phenomenological couplings as 
functions of /? we have computed the Taylor series of these quantities around Ps up 
to the third order, where in the simulation the configurations are generated with a 
Boltzmann factor corresponding to jSg- Since the Taylor expansion provides accurate 
results only in a small neighbourhood of Ps, we tried to chose Ps ~ Pkt{Li), where 
Pkt{.Li) is the solution of eqs. fl32l) for the given lattice size Li. To this end, we 
performed some preliminary simulations or we extrapolated the result for Pkt{Li) 
from small lattice sizes Li to large ones. We carefully checked that in our final 
simulations the differences Pkt{Li) — Ps are small enough to ensure that the Taylor 
series evaluated at Pkt{Li) are sufficiently accurate for our purpose. 

In our final simulations, mostly, we performed 2 x 10^ measurements. Only in 
the case our largest lattice, 32 x 1024^, we only measured 1.4 x 10^ times. For 
each measurement, we performed single cluster updates, boundary fiip updates and 
3 overrelaxation sweeps and one Metropolis sweep. The number of single cluster 
updates was chosen such that, on average, the lattice is covered once by the clusters. 

E.g. for the 32 x 1024^ lattice, the integrated autocorrelation time of the mag- 
netic susceptibility is ~ 0.7. I.e. individual measurements are almost statistically 
independent. The total CPU time used for these simulations was roughly the equiv- 
alent of 9 years on a single core of a 2.6 GHz Opteron CPU. Our results obtained 
from the matching with the dual of the ASOS model at /? = 0.80608 are summarized 
in table [5l 

For all values of Lq, Pkt{Li) is decreasing with increasing Li. Apparently 
the results are converging; The difference Pkt{Li) — Pkt{'2Li) is decreasing with 
increasing Li. For Lq = 4 and 6, where we have the largest range of lattice sizes 
available, the differences of the results for Li/Lq = 64 and Li/Lq = 128 are already 
smaller than the error-bars. 

Eqs. (1321) are supposed to give the universal behaviour within the numerical 
precision that has been reached in ref. [M]. Corrections to the universal behaviour 
are expected to decay with a power of Li, possibly multiplied by some power of 
InLi. There are certainly corrections with an exponent e = 7/4 due to the analytic 
background in the magnetic susceptibility. Furthermore there are corrections due 
to irrelevant scaling fields. 

In order to check these expectations, we have fitted our estimates of Pkt{Li) for 
Lq = 4 and 6, where we have the largest values of Li/Lq available, with the ansatz 

Pkt{L,) = fc(oo) + cL^' , (34) 

where Pkt{oo), c and e are the free parameters of the fit. The results of these fits 
are summarized in table [61 The numerical results for the correction exponent e are 
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Table 5: Matching the thin film with the dual of the ASOS model at ^ = 0.80608, 
eqs. ( f32j) . The parameters of the matching are the inverse KT temperature 13 kt of 
the thin film and the lattice size L^sos of the dual of the ASOS model. 
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192 


0.5321416(48) 


65.9(6) 


12 


384 


0.5321010(33) 


139.6(1.8) 


12 


768 


0.5320864(28) 


297.6(6.6) 


16 


128 


0.5245609(34) 


33.87(13) 


16 


256 


0.5244956(27) 


67.5(5) 


16 


512 


0.5244626(23) 


140.5(1.8) 


16 


1024 


0.5244547(19) 


296.9(6.1) 


24 


192 


0.5177878(19) 


34.39(13) 


24 


384 


0.5177446(15) 


68.48(48) 


24 


768 


0.5177340(11) 


143.6(1.6) 


32 


256 


0.5148522(12) 


34.88(13) 


32 


512 


0.5148259(9) 


70.4(4) 


32 


1024 


0.5148147(9) 


147.7(2.2) 
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Table 6: Results of fits with the ansatz ([31]). Li^rnin is the smallest lattice size that 
has been included into the fit. For a discussion see the text. 





f-'l,'min 


(3kt{oo) 


c 


e 


xVd.o.f. 


4 


32 


0.609669(11) 


0.14(4) 


1.47(9) 


2.78 


4 


64 


0.609681(13) 


0.6(6) 


1.8(3) 


4.31 


6 


48 


0.568241(8) 


0.12(4) 


1.40(9) 


1.48 


6 


96 


0.568250(8) 


0.6(8) 


1.8(3) 


0.94 



Table 7: Matching the thin film with the dual of the ASOS model at (3 = 0.80606, 
eqs. ( 133)) . As examples we give the results for Lq = 4 and 6. 





Li 




Lasos 


4 


32 


0.610532(24) 


28.08(9) 


4 


64 


0.609996(18) 


54.1(3) 


4 


128 


0.609793(17) 


114.0(1.3) 


4 


256 


0.609695(12) 


227.7(3.2) 


4 


512 


0.609705(11) 


467.(12.) 


6 


48 


0.568780(14) 


30.57(10) 


6 


96 


0.568459(12) 


60.59(36) 


6 


192 


0.568318(8) 


125.3(1.4) 


6 


384 


0.568266(7) 


260.1(4.6) 


6 


768 


0.568262(7) 


542.(18.) 



not very precise; they are essentially consistent with a leading correction exponent 
e = 7/4 due to the analytic background in the magnetic susceptibility. 

In order to check the uncertainty of our results for Pkt caused by the uncertainty 
of the estimate of the KT transition temperature in the ASOS model, we have 
repeated the matching with eqs. ( l33l) . As examples we give the results obtain for 
Lq = 4 and Lq = 6 in table [71 

We find that the difference in Pkt compared with the matching with eqs. (132!) 
is smaller than the statistical error of Pkt- Also in the case of the matching lattice 
size Lasos of the ASOS model the difference is small. These observations also hold 
for larger values of Lq not given in table [71 
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Table 8: Final results for the inverse KT transition temperature I3kt of thin films 
of the thickness Lq. The statistical error is given in () while the systematic one is 
quoted in []. For a discussion see the text. In addition we give in the last row the 
inverse critical temperature of the three dimensional system. 





f^KT 


4 


0.60968(1)[1] 


6 


0.56825(1)[1] 


8 


0.549278(5) [9] 


12 


0.532082(3) [5] 


16 


0.524450(2) [3] 


24 


0.517730(2) [2] 


32 


0.514810(1)[2] 


oo 


0.5091503(6) 



As our final estimate of the inverse KT transition temperature for Lq = 4 we 
take Pkt = 0.60968(1) [1], which is the result of the fit with the ansatz using 
Li,Tnin = 64. The number given in [] is an estimate of the systematic error. It is 
estimated from the difference of the fit with the ansatz (|34l) using Li^rnin = 32 and 
Li,min = 64. Furthermore we take into account the difference between the matching 
with eqs. ( I32i) and eqs. (l33l) . In a similar way we arrive at Pkt = 0.56825(1) [1] for 
Lo = 6. 

For Lo > 8 the range of Li/ Lq that we have simulated is smaller than for Lq — 4 
and 6. Therefore, for these values of Lq, we use e = 7/4 fixed in our fits ( l34l) to 
obtain (3kt in the limit Li/Lq oo. In all cases we did not include the smallest 
lattice size, Li/Lq = 8, into the fit. Even if e = 7/4 is the dominant correction, 
subleading corrections will cause systematic errors. To estimate these, we assume 
that the effective correction exponent might take a value in the range 1.5 < e < 2 
as it is suggest by our fits for Lq = 4 and 6. To this end, fitting with the ansatz (!34l) 
we have used e = 1.5 in addition to e = 1.75. The systematic error is then estimated 
by the difference of these two fits. As a further check of systematic errors, we have 
repeated these fits for the estimates of Pkt obtained from the matching with the 
dual of the ASOS model at /3 = 0.80606, eqs. fl33|) . Our final results are summarized 
in table E 

From the matching of the thin films with the dual of the ASOS model, we obtain 
also a scale factor that relates the lattice size of the thin films and the dual of the 
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Figure 3: In the figure we give the matching factor (Lq + Ls)Lasos/Li between 
thin films and the dual of the ASOS model. 

ASOS model. Following the RG prediction, thin films of different thickness Lq 
match, up to corrections to scaling, for the same values of Li/{Lq + Lg). 
In figure [3] we have plotted the matching factor 

bAsosjiim = Lasos{Lo + Lg)/ Li (35) 

between the dual of the ASOS model and thin films, where we have used = 1.02, 
eq. ( !26|) . Starting from Lq = 6 the results for the same Li/Lq are consistent among 
different Lq, confirming the RG prediction. The results for the scale factor are all 
in the range 4.2 < hAsosjUm < 5.2. Unfortunately, the statistical error increases 
a lot with increasing Li/Lq. Also the value of the scale factor bAsosjUm is still 
increasing for our largest values of Li/Lq. Taking into account this trend, we shall 
assume 

bAsosjiim = 5.0(5) (36) 
as estimate for the limit Li/Lq oo in the following. 
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Using the matching factors between the dual of the ASOS model and other XY 
models that are given in ref. |19], we arrive at 



bxYJilm = LxyLqI L 
bvillainjilm = LviUainLo/ L 

bBcsosjilm = LbcsosLq/ L 



1 



1 



1 



1.7(2) 

0.58(6) 

1.8(2) 



(37) 



In the following we like to further check that the thin films indeed undergo a 
Kosterlitz-Thouless transition. In figurelHwe give [/4 and ^2nd/L at the estimates of 
(5kt given in table[Hlas a function of Li/ (Lq + Ls), where we use Ls = 1.02, eq. ( l26l) . 
For comparison we also give the results for the dual of the ASOS model represented 
by eqs. (132|) . In this case, f/4 and ^2nd/ L are plotted versus LAsos/bAsosjUm, using 
bAsosjiim = 5. The data for thin films with different Lq fall nicely on top of each 
other. This is even true for our smallest values of Li/Lq. The thin film results do 
match nicely with those of the ASOS model only for the larger values of Li/Lq. 
The approach of the thin film results to those of the ASOS model seems compatible 
with power law corrections. This behaviour is, of cause, directly related with the 
behaviour of (3kt{.Li) obtained from the matching of t/4 and i2nd/L. Also there 
we have seen sizable corrections for small values of Li/Lq that could be fitted by a 
power law. 

In figure [5] we give our results for the ratio of partition functions Za/Zp and 
the helicity modulus T. In contrast to ^2nd/ L and f/4 these quantities have not 
been used to determine [3kt- Therefore they provide an additional check of the 
rehability of our results for [3kt and the matching factor bAsosjUm- Also in the 
case of Za/Zp the data for different values of Lq fall reasonably well on top of each 
other. Unfortunately we have no data for Za/Zp for a two-dimensional XY model 
available. Therefore we can only compare with the leading behaviour, eq. fl69|l . at 
the KT transition that we derive in the appendix. We determined the constant C 
by a 1-parameter fit of our data for Lq = 16 using Li = 512 and Li = 1024, where 
we have used L = Li/{Lq + Lg). The result is C = 4.31(7). The result of this fit is 
given in the upper part of figure O by the dashed purple line. 

In the lower part of figure [5] we give L^T . Note that, in contrast to section \T?2\ 
the T used here is measured on the thin films. The results for different Lq fall nicely 
on top of each other. As expected, the relative error of LqT increases with increasing 
thickness Lq of the film. For comparison we also plot results obtained for the 2D XY 
model [311 . We have plotted results for Lxy = 12, 16, 24, 32, 48, 64, 96, 128, 192 and 
256 as a function of L = Lxy/bxYjum- The dashed hne only connects the points 
to guide the eye. Here we find an almost perfect matching between the thin films 
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gure 4: Effectively two-dimensional QAite size scaling of f/4 (upper part of the 
;ure) and C,2nd/L (lower part of the figure). For a discussion see the text. 



and the two-dimensional XY model, starting from our smallest values of Li/Lq. 
This suggests that corrections to the universal KT behaviour are very small in the 
helicity modulus. 



7.1 Scaling of (3kt with Lq 

As the thickness Lq increases, the KT transition temperature of the thin film ap- 
proaches the transition temperature of the three-dimensional system. Finite size 
scaling predicts that this approach is governed by eq. ([3]), or equivalently in terms 
of the inverse temperature 

(5kt{L^) - Pc,3d ^ L-^"" . (38) 

Indeed from table [8] we read off that j3KT is approaching /3c,3d as Lq is increasing. 
In order to check whether this approach is compatible with a power law and in 
particular with the predicted exponent, we have computed 

^effALo) = - ln(2)/ In - — -— , (39) 

V PKT{-I^0) — Pc,3d / 

where we use Pc,3d = 0.5091503(6) |T0]. Analogously we define h'eff,T{LQ), where in 
eq. fl39|) Pkt — Pc,3d is replaced by TKT — Tc^3d- In the limit Lq — »• cxo both definitions 
give, by construction, the same result. Our results based on the numerical estimates 
of (3kt given in table [H] are plotted in figure El Note that the error-bars are smaller 
than the symbols. For simplicity we have added the systematic and the statistical 
error that is given in table [H Also the error due to the uncertainty of /?c,3d is taken 
into account. 

Our results for u^ffA^o) and Peff^riLo) differ quite a lot, indicating that analytic 
corrections are still important for the values of Lq considered. In both cases the 
effective exponent is decreasing with increasing Lq. Even for the largest value of 
Lq, the effective exponent is by about 4% and 6% larger than the asymptotically 
expected v = 0.6717(1). 

In the case of the improved model, leading corrections to scaling are due to 
the free boundary conditions. To take these into account we define an improved 
effective exponent as 

f 2L„ + L, \ f !3kt(2L„) - 3,,m \ 



Also here we define z/e//,T,L,(-^^o) by replacing (3kt - Pc,3d by Tkt - T^m in eq. 
Using Ls = 1.02, eq. fl26|) . we arrive at the results that are plotted in figure |6l 
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Table 9: Results of fits with the ansatz ([HP. Lo^min is the minimal thickness that 
is taken into account. As discussed in the text, a, Lg and c are the parameters of 
the fit. 





a 


Ls 


c 


xVd.o.f. 


4 


1.0321(11) 


1.158(29) 


1.38(11) 


3.80 


6 


1.0286(13) 


1.037(36) 


0.90(13) 


0.39 


8 


1.0284(18) 


1.030(59) 


0.87(23) 


0.57 



Table 10: Similar to table\B but with the ansatz 



f-'0,min 


a 


Ls 


c 


xVd.o.f. 


4 


3.9788(25) 


1.118(11) 


-0.851(24) 


2.43 


6 


3.9697(35) 


1.057(21) 


-1.03(6) 


0.43 


8 


3.9684(55) 


1.047(37) 


-1.07(11) 


0.60 



Now, i'eff,i3,Ls is increasing with increasing Lq, while i'eff,T,Ls is decreasing. Now, 
for the largest value of Lq, the difference with u = 0.6717(1) is only about 1% in 
both cases. The remaining difference between i'eff,/3,Ls and i^eff,T,Ls indicates that 
analytic corrections have to be taken into account. Therefore we have fitted our 
data for /Skt with the ansatz 

fc(^o) - /3c,3D = a{Lo + Ls)-^/'' X (1 + c{L, + 1^)-^/'') , (41) 

where we have included the leading analytic correction in addition to the correction 
caused by the boundaries. We take a, Ls and c as free parameters of the fit, while 
now /3c,3D = 0.5091503 and v = 0.6717 are fixed. I.e. here we assume the RG 
prediction for the exponent to be correct. We have not included corrections oc L^'^ 
since for A = 2.1 the amplitude of these corrections should be small. Furthermore 
it should be difficult to disentangle them from the boundary corrections which are 
oc Lq^. The results of these fits are given in table O Next we have repeated the fit, 
replacing /3 by T in the ansatz fKT]) : 

Tc,3D - Tkt{Lo) = a{Lo + LsY^/" x (1 + c(Lo + LsY^'") . (42) 

For both types of fits we find x^/d-o.f. < 1 already for Lo^min = 6, where we 
have included all data with Lq > Lo^min into the fit. The result for Ls from the fits 
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with Lo^rnin = 6 and 8 are compatible among each other. Furthermore they are in 
perfect agreement with the result = 1.02(7), eq. fl26l) . obtained from FSS at Pc,3d- 
Also the results for the leading amplitude are compatible among the two types of 
fits: Using a = a(3l we convert the results of a given in table [10 to a = 1.0314(6), 
1.0291(9) and 1.0287(14) for Lo^min = 4, 6 and 8, respectively. I.e. the difference 
between the two types of fits is negligible. 

We take our final result for the amplitude from the fit with ansatz fHT!) and 
Lo^min = 8. We have repeated this fit for shifted values of the input parameters: 
for /3c,3d = 0.5091509, z/ = 0.6717 and for (5^,^ = 0.5091503, u = 0.6718. From the 
results of these fits we conclude for the amplitude 

a = 1.028(2) - 10.5 x (z/ - 0.6717) - 400 x (/J^.sd - 0.5091503) . (43) 

The analysis of the behaviour of (3kt presented in this section leaves little doubt 
that the finite size scaling prediction (13|38!) is correct. To fit data for thicknesses 
in the range 6 < Lq < 32 properly, boundary corrections and leading analytic 
corrections have to be taken into account. Note that in the case of models which 
are not improved, also corrections oc L^'^ should play a role. 

7.2 Universal amplitude ratios 

As we have noted already in the introduction, finite size scaling predicts that the 
KT transition should occur at a universal value of Lq/^, where ^ is a correlation 
length of the three-dimensional bulk system. In this section we shall compute this 
ratio, using the transversal correlation length = 1/T and the second moment 
correlation length ^2nd in the high temperature phase. 

For the following discussion it is more convenient to consider the inverse of the 
function Pkt{Lq), namely the thickness Lq^kt{P) where the KT transition occurs 
at a given value of /?. Of course this is not well defined, since Lq can take only 
integer values, and therefore for almost all values of f3 there is no such Lq. However 
by proper interpolation this problem can be solved, at least on a practical level. 

There are several strategies to compute the universal amplitude ratios from 
Monte Carlo data. For a discussion see e.g. ref. [50]. If both quantities live in 
the same phase, as it is the case here for Lq^kt and T, we can simply consider the 
product at the same value of /?, or equivalently the same temperature T. Then the 
critical limit is taken 

[Lo,KTrY= lim Lo,KT{f3)T{(3) . (44) 
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One expects from RG theory that the approach to the critical hmit is described by 
(Lo,^T + L,)T = [Lo,KTrr + cL.^^KT + - ■ (45) 

If the two quantities hve in different phases, as it is the case here for Lq^kt and C,2nd 
in the high temperature phase, we define 

[Lo,KT/^2ndT = Jim Lo,i,T(/3)/6nd(2/3c,3d - /5) • (46) 

I.e. the two quantities are taken at the same distance from the critical point of 
the three-dimensional system. This is however not uniquely defined. Instead of the 
same distance in j3 we could also use the same distance in T. Therefore, in contrast 
to eq. ( I45l) also analytic corrections appear 

The universal ratio can also be expressed in terms of amplitudes: 

[Lo,KT/^2ndr = a^f2nd,+ , (48) 

where a is defined by eq. pil) and f2nd,+ is the amplitude of the second moment 
correlation length in the high temperature phase. 

7.2.1 The second moment correlation length in the high temperature 
phase 

Here we analyse the Monte Carlo data for ^2nd given in table 5 of ref. fT2]- We have 
fitted these data with the ansatz 

i2nd = f2ndA(3c,3d ' f^Y'' X [1 + 6(/5,,3rf - (3^ + c(/5,,3d - P)] , (49) 

where we have included leading corrections oc (/5c — /?)^ and analytic corrections. 
We have fixed v = 0.6717, p = 0.5091503 and A = ucu = 0.527, given in ref. 
The results of these fits are summarized in table [TTJ A x^/d.o.f. close to one is 
reached for /3mm = 0.485. The result for f2nd,+ shows very little dependence on 
Pmin- Therefore we take as our final result the one for Pmin = 0.49: 

f2nd,+ = 0.26362(8) + 223 x (/5c,3d - 0.5091503) - 2.1 x (z/ - 0.6717) , (50) 

where we have obtained the dependence on /3c,3d and on u by redoing the fits with 
slightly changed values of Pc,3d and z/. 
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Table 11: Fitting ^2nd in the high temperature phase of the three-dimensional 0'^ 
model at X = 2.1 with the ansatz ( BPj) . All data with {3 > fimin Sire taken into 
account. 





f2nd,+ 


b 


c 


xVd.o.f. 


0.48 


0.26352(5) 


0.050(4) 


-0.77(2) 


1.54 


0.485 


0.26359(6) 


0.043(5) 


-0.73(2) 


1.14 


0.49 


0.26362(8) 


0.039(8) 


-0.72(4) 


1.23 



Now we can compute the universal ratio by using the results fl43ll50p : 
[Lo,KT/^2ndT = «'/2;d,+ = 3.864(6) -4300(/?e,3d- 0.5091503) + 5(z/ - 0.6717) . (51) 

Note that the error of this product is mainly due to the error of a. Taking into 
account the uncertainty in /3c and u we arrive at 

[i^o,W6nd]* = 3.864(9) . (52) 

One should note that the exponential correlation length C,exp, which describes the 
asymptotic decay of the correlation function, differs only slightly from C,2nd which 
has been used here. Following ref. [lOj: 

lini 1^ = 1.000204(3) , (/3 < /3c,3d) ■ (53) 

Hence at the level of our accuracy this difference plays no role. 

7.2.2 The helicity modulus in the low temperature phase 

The results for the helicity modulus, which are summarized in the second column 
of table [121 are taken from ref. [T2]. In order to obtain Lq^kt{i3) at the values of (3 
that were simulated in ref. [12] , we have to interpolate the results given in table [81 
To this end we have used the ansatz (HTl) along with the values for the coefficients 
obtained from the fit with Lornin = 8. The numbers for Lq xt{P) that are given 
in the third column of table [T^l are obtained by numerically inverting eq. (I4ip . 
Note that for the present purpose, the asymptotic correctness of the coefficients of 
eq. fHTl) is less important. For the present purpose it is sufficient that the function 
accurately describes our data in the range of Lq that has been simulated. In table 
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Table 12: Results for the product Lo^kt{P)'^{P) 



p 


T 


Lo,KT 




{Lqkt + Ls)T 


0.515 


0.04939(8) 


31.29 


1.5454(25) [8] 


1.5958(26)[8]{35} 


0.52 


0.07456(7) 


20.37 


1.5187(14)[8] 


1.5948(15)[8]{52} 


0.525 


0.09628(7) 


15.61 


1.5029(11)[8] 


1.6011(12)[8]{67} 


0.53 


0.11565(7) 


12.845 


1.4855(9) [7] 


1.6035(10)[8]{81} 


0.535 


0.13349(6) 


11.012 


1.4670(7) [7] 


1.6062(7) [8]{93} 


0.54 


0.15035(6) 


9.692 


1.4572(6) [7] 


1.6105(6)[8]{105} 


0.55 


0.18145(6) 


7.895 


1.4325(5) [7] 


1.6176(5)[8]{127} 


0.58 


0.26340(5) 


5.229 


1.3773(3) [7] 


1.6460(3)[8]{184} 



fT2]we quote no error for Lq^kt- Given the accuracy of Pkt{Lq) we expect an relative 
error of about 1/2 %o for Lq^kt at a fixed value of [3. 

In the fourth column of table [12] we have computed the product Lq^kt{i3)'^{i3)] 
in we give the error due to the error of T and in [] the 1/2 %o error assumed 
for Lo^kt{P)- -Z^o,A'r(/3)T(/?) is increasing with decreasing /?. Even for the two 
smallest values of the inverse temperature /3 = 0.52 and /3 = 0.515, the difference of 
Lo,kt{P)'^{P) is larger than the error. In order to take into account the boundary 
corrections, which are leading for the improved model, we have computed {Lq kt + 
Ls)T. The results are given in the last column of table [121 We have used the 
numerical value Lg = 1.02(7), eq. (1261) . The error of Lg is refiected by the number 
given in {}. Now the dependence of the product on /3 is much reduced. Starting 
from j3 = 0.53, the results are consistent within the quoted errors. As our final 
result we take 

[Lo,KTrY = 1.595(7) (54) 

which is the average of the result for (3 = 0.515 and 0.52. The error-bar includes all 
results up to /? = 0.53. For illustration we have plotted in fig. [7] the results given 
table [12] along with our final result (JSl). 



32 



i.6.r 



I i- 



1 i i 



i4 1.5 h 

I 



1.4 - 



L_0,KT Y 
(L_0 + L_s) Y 



1-3, 



0.05 



0.1 
1/L_0,KT 



0.15 



0.2 



Figure 7: The numerical estimates of Lq^kt'^ and [Lq^kt + Ls)T which are quoted 
in table [12] are plotted as a function of In addition we give our final 

result (|54|) (green diamond) for the limit Lq^kt ^ oo. For a discussion see the text. 
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Table 13: Preliminary estimates of (3kt for the dynamically diluted XY model at 
D = 1.02. The results are all obtained from lattices of the size Li = L2 = 16 x Lq 
by requiring ^2nd/ L = 0.792(2). For a discussion see the text. 







4 


0.6767(3) [2] 


8 


0.60871(14) [6] 


16 


0.58083(11) [3] 


32 


0.57012(3) [2] 



8 Preliminary results for the dynamically diluted 
XY model and the standard XY model 

We have performed a preliminary study of the dynamically diluted XY model at 
D = 1.02 and the standard XY model. To this end we have simulated lattices of 
the size Li = L2 = 16 x Lo with Lq = 4, 8, 16 and 32. In order to determine Pkt 
we assume that effectively two-dimensional finite size effects are universal in thin 
films in the 3D XY universality class. I.e. that they are the same as for the 0^ 
model studied in section [71 In particular we read off from figure S] that ^2nd/L = 
0.792(2) for Li/Lq = 16 at the KT transition temperature. In the following we 
have determined an estimate of Pkt by requiring that (2nd/ L assumes the value 
0.792 on lattices with Li/Lq = 16. The results obtained this way are summarized 
in table US] for the dynamically diluted XY model at D = 1.02 and in table for 
the standard XY model. The number given in () is the statistical error, while the 
number given [] is the error due to the uncertainty of ^2nd/L. 

As a consistency check we have computed LqT and U4 at the values of Pkt 
summarized in tables [T3] and UM Indeed we find that U4 ~ 1.02 and LqT ^ 0.715 as 
we also read off for Li/Lq = 16 from figures H] and [5] for the 0^ model, respectively. 

The results for the ddXY model, including all data, can be fitted by 

ddXY,D = 1.02 : pKT = 0.5637963 + 1.1286(71) x (Lo + 0.70(3))"^/°-*^^^^ . (55) 

Since we have only data for a few thicknesses of the film, we cannot check for 
systematic errors due to subleading corrections. 

In table [HI for the standard XY model we also give the results that can be found 
in table I of ref. [30] which were obtained with staggered boundary conditions. 
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Table 14: Preliminary results for (3kt for the standard XY model are summarized in 
the second column. In the third column we give for comparison the results obtained 
in ref. J^D|/ . Note that these refer to staggered boundary conditions and therefore 
no exact match of the numbers is expected. 





Pkt, ours 


pkt, m 


4 


0.5665(3) [3] 


0.5448(3) 


8 


0.4987(2) [1] 


0.4949(10) 


12 




0.47934(23) 


16 


0.47132(3) [3] 


0.47226(7) 


20 




0.46847(15) 


32 


0.46050(2) [1] 





Therefore no exact match with our results is expected. For Lq = 4 the inverse KT 
transition temperature i3kt of ref. [30j is clearly smaller than ours. For Lq = 8 it 
is still slightly smaller, while for Lq = 16 it is slightly larger than ours. 

Fitting the results for the standard XY model is more difficult than those of the 
ddXY model at D = 1.02, since for the standard XY model we expect a sizable 
correction oc with u ~ 0.785(20) [in]. Ignoring these corrections we get from 
fitting all our data of table [14] 

XY : (3kt = 0.4541655 + 1.1387(45) x (Lq + 0.75(2))-^/°-^^^^ (56) 

with x^/d.o.f. = 2.02. Taking into account the correction oc L"*^ instead of the 
boundary correction we obtain 

XY : (3kt = 0.4541655 + 1.1545(46) x (1 - 0.70(2) x L^^-^^^) x (57) 

with x^/d.o.f. = 1.92. Given the small number of data points available we made 
no attempt to fit both types of corrections simultaneously. The difference between 
the two fits fl56f57l) might give some indication of systematic errors. 
Fitting the values of Pkt for Lq = 8, 12, 16 and 20 of ref. [30] we get 

PKTMaggered = 0.4541655 + 1.47(13) X (Lq + 3.13(76))"^/°-^^^^ (58) 

or alternatively, taking into account corrections oc Lq^: 

PKT,staggered = 0.4541655 + 1.44(10) X (1 - 1.92(25) X V'^^^) X L,'^'''''' . (59) 
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If free and staggered boundary conditions are equivalent, the prefactor of Lg 
should be identical. Apparently, here we find a discrepancy of about three times 
the error. Based on the data of ref. [30] it is hard to decide whether this discrep- 
ancy is due to subleading corrections which are not taken into account or a subtle 
cancellation of corrections oc Lq*^ and oc Lq ^ or that free and staggered boundary 
conditions gainst our expectation, not equivalent. 

In relation with ref. [TO] we have also determined the correlation length of the 
standard XY model and the ddXY at D = 1.02 in the high temperature phase. 
Using a similar analysis as for the 0^ model at A = 2.1 in the previous section, we 
arrive at /2„d,+ = 0.2822(3) for the ddXY model at D = 1.02 and /2„d,+ = 0.2880(3) 
for the standard XY model, where the amplitude f2nd,+ is defined by eq. ( H^ . With 
eq. ( I5T]) we arrive at [Lq^kt / ^2nd]* = 3.84(2) using eq. (155]) for the ddXY model at 
D = 1.02 and [Lo,KT/^2ndY = 3.79(1) using eq. m and [Lo,KT/6nd]* = 3.82(1) 
using eq. fl57|l for the standard XY model. Given the fact that systematic errors are 
not fully taken into account, we regard these results as consistent with our more 
precise estimate (1521) obtained from the (p^ model at A = 2.1, confirming universality. 



9 Comparison with experimental results 

To our knowledge, there are no previous theoretical results for the amplitude ratios 
[Lo^KT^Y and [Lo,i^T/6nd]*- 

The experimental study of ^He provides no access to the correlation length in the 
high temperature phase. On the other hand the helicity modulus can be computed 
from the superfluid density. Following ref. [12] the helicity modulus as defined here 
is given by 

'^(T) = , (60) 

where PsbiT) is the superfluid density of the bulk system and m the mass of a "^He 
atom. In the literature the inverse of the helicity modulus is called transverse corre- 
lation length = 1/T. The superfluid density can be obtained from measurements 
of the specific heat and the velocity of the second sound. 

A number of experimental studies of the KT transition in thin films of ^He have 
been performed. For an overview see ref. [201 • Only a few of these works give 
results for the KT temperature, obtained e.g. from the onset of the superfluidity, 
for a sufficiently large range of the thickness of the film to extract a result for 
the universal combination [Lq^a't'^]*- In fact in two of them an explicit result for 
[Lq^kt'^Y is quoted: 
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Sabisky and Anderson [51] have performed experiments with films of ^He on 
various substrates. They have studied films of a thickness up to about 75 A. For 
films on CaF2 they find for Lq ~ 35 A power law scaling of the film thickness 
with the reduced KT temperature with an exponent 2/3, which is in reasonable 
agreement with the correlation length exponent v = 0.6717(1) [lOj. Concerning the 
scaling of the thickness of the film with the superfiuid density of the bulk system 
they write "the product {ps/p)C,(T) = 5.8T over the entire temperature range". In 
their notation ^(T) is the thickness of the film where the KT transition takes place 
at the temperature T. It is clear from the text that Sabisky and Anderson have 
omitted a factor A/K on the right hand side of the equation. Plugging in numerical 
values for p, h, m and Tx into eq. (!60|) we arrive at 

[Lq^ktT]* = 5.8 X 0.2664 = 1.545 . (61) 

The authors of ref. |52j have studied Helium films on glass of a thickness up 
to 47 A. They also find power law scaling of the KT thickness of the film with 
the temperature. Their numerical estimate for the exponent is 0.71. Even more 
interesting, they quote in their eq. (7): 

4(T) = 0.43 + 1.6ia(T) , (62) 

where Sc{T) is the KT thickness of the film at the temperature T and ^±(T) is the 
transverse correlation length of the bulk. I.e. they give a direct estimate for the 
universal combination which we are interested in. It matches perfectly with our 
theoretical result. 

It is a bit puzzling that there are more recent experimental studies using larger 
thicknesses of the ^He film, which fit less well the theoretical expectation than the 
two reported above. Yu et al. [53] have studied ^He films on Kapton up to a 
thickness of 156 A. They see a power law scaling, starting from 25 A. However they 
find an exponent 0.52(1) which is clearly smaller than u = 0.6717(1). Furthermore, 
from their figure 3 it can be clearly seen that the ratio of ^±(T) and the KT thickness 
varies a lot over the range of thicknesses that has been studied. For their largest 
thickness, 156 A we read off from their figure Lq^kt/^± ~ lO'^'^ ~ 1.58. On the 
other hand, for 25 A one reads off Lq^kt/^± ~ 10°'^ ~ 2.51. The reduced transition 
temperatures are 0.2541 and 0.0073 for 25 A and 156 A , respectively |54j . 

In refs. [5^ the KT temperature has been determined for films of the thick- 
ness 483 A and 2113 A at saturated vapour pressure. These films are confined 
by two wafers of silicon. The adiabatic fountain resonance method has been used 
to determine the KT transition temperature. The reduced transition temperatures 
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are ^ 1.29 x lO^^ and ^ 1.64 x 10"^ for 483 A and 2113 A , respectively [M]. 
An estimate of can be obtained e.g. from ref. [7j: from their figure 3 we 
read off .^o ~ 3.42 A at saturated vapour pressure, where ~ ^o^~^- It follows 
Lo,kt/U ~ 1-62 and 1.77 for 483 A and 2113 A , respectively. 

10 Summary and Conclusions 

We have simulated lattice models in the three-dimensional XY universality class 
with thin film geometry. In order to mimic the vanishing order parameter at the 
boundaries that is observed in experiments on thin films of ^He we have implemented 
free boundary conditions. These boundary conditions lead to corrections to finite 
size scaling cx: Lq^, where Lq is the thickness of the film. These corrections can 
be described by Lo^^ff = Lq + L^. Furthermore one expects corrections to finite 
size scaling oc Lq'^ [TH], where u = 0.785(20) [lU]. Analysing Monte Carlo data, it 
is difficult to disentangle corrections with exponents that are so close. In fact in 
previous studies, e.g. [30], of thin films, using the standard XY model, only the 
boundary corrections have been taken into account, which might lead to systematic 
errors in the extrapolation to the critical limit. In order to avoid this problem, we 
have studied improved models. In particular we have studied in detail the 0^ model 
on a simple cubic lattice at A = 2.1. This value of the parameter A is close to the 
most recent estimate A* = 2.15(5) [10], where A* is defined such that the amplitude 
of leading corrections exactly vanishes. At A = 2.1 the amplitude of corrections to 
scaling should be at least by a factor of 20 smaller than in the standard XY model. 

First we have studied the 0^ model at the critical point of three-dimensional 
system. We have simulated lattices of the size Lq = Li = L2 and Lq + 1 = Li = L2 
with Li = 8, 12, 16, 24, 32 and 48. We employed free boundary conditions in 0- 
direction, while periodic boundary conditions are employed in 1 and 2-direction. 
We have measured the four phenomenological couplings U4, Uq, ^2nd/L and Za/Zp. 
The latter two quantities are taken only for the 1 and the 2-direction. We clearly 
see corrections oc L^^ in all of these quantities. This is in clear contrast with 
simulations, where periodic boundary conditions in all directions are employed. E.g. 
in ref. using improved models one sees corrections oc Lq^, where 1.6 ^ e ^ 2. 
Casting the corrections oc Lq ^ into the form Lo,e// = Lq + Lg, one gets, within the 
error, the same value for Lg from all four quantities. As our final result we quote 
Ls = 1.02(7). Note that depends on the model and on the particular boundary 
conditions that are used. In future studies of other properties of thin films, like 
e.g. the thermodynamic Casimir force, using the same lattice model as here, the 
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knowledge of Lg is a valuable asset. 

Next we studied the phase transition in thin films. To this end, we applied a 
finite size scaling method that we have proposed recently [34] . Finite size scaling 
of the effectively two-dimensional film means Lq <C Li,L2 ~ ^, where ^ is the 
correlation length in the 1 or 2-direction of a system with the thickness Lq and 
infinite extension in the other two directions. Also here, we have studied in detail 
the (f)^ model at A = 2.1. We have simulated lattices up to Li/Lq = 128 for 
Lq = 4 and 6, up to Li/Lq = 64 for Lq = 8, 12 and 16 and up to Li/Lq = 32 
for Lq = 24 and 32. Matching the finite size behaviour of the second moment 
correlation length C,2nd/L and the Binder cumulant U4 with that of the dual of 
the ASOS model at the Kosterlitz-Thouless (KT) transition we obtain accurate 
estimates of the transition temperature of the thin films. Note that the finite size 
scaling behaviour is very similar for different two-dimensional XY models [M] . The 
dual of the ASOS model has been chosen in ref. [31] as reference model for purely 
technical reason. Furthermore, we obtain an estimate for the scale factor hAsosjUm 
that relates the effective two-dimensional lattice size Li/Lq of the thin film and that 
of the dual of the ASOS model. The finite size scaling behaviour of Za/Zp and the 
helicity modulus times the thickness LqT, which were not used to determine I3kt 
and bAsosjiim, clearly confirm the KT nature of the transition in the thin films. 

Next we have studied the scaling of the inverse of the KT temperature Pkt 
with the thickness Lq of the film. We find that for the thicknesses that we have 
studied, even in the improved model, corrections to scaling are important. Using 
the ansatz (HT!) that includes corrections due to the free boundary conditions as 
well as leading analytic corrections we can fit, fixing u = 0.6717(1), all data with 
Lo > 6, with xVd.o.f.=0.39. 

For the first time, we provide a theoretical estimate for the amplitude ratio 

[Lo,KT/^±r = 1-595(7) , (63) 

where C,± is the transversal correlation length in the low temperature phase of the 
three-dimensional bulk system. Note that C,± is the inverse of the helicity modulus 
as defined here. This result nicely compares with the experimental results of refs. 
[51] [52] . On the other hand, more recent experiments [55l [56i| on ^He films with 
thicker films seem to match less well. 

We have also computed the universal amplitude ratio [Lq^kt / ^2nd]* = 3.864(9), 
where C,2nd is the second moment correlation length of the three-dimensional bulk 
system in the high temperature phase. We also provide some preliminary results for 
the KT transition temperature for the dynamically diluted XY model at D = 1.02 
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and the standard XY model for Lq = 4, 8, 16 and 32. Using the data of these two 
models wc get results for [I/o,KT/^2nd]* that are essentially consistent with the one 
quoted above for the 0^ model, confirming the universality of the result. 



11 Appendix A: Ratio of partition functions ZajZ^ 
at the KT transition 

In the spin-wave approximation of a two-dimensional XY model, the partition func- 
tion is given by 

Zsw= J] Vr(ni,n2)Z(0,0) , (64) 

ni,n2 

where n\ and n-i count the windings of the XY field along the 1 and 2 directions. 
In the case L\ — L2, for isotropic couplings, the weights are given by 

W{ni,n2) =exp{-27i^Psw[nl + nl]) , (65) 

where for periodic boundary conditions, ni and n2 take integer values and in the 
case of antiperiodic boundary conditions half-integer. I.e. in the spin-wave approx- 
imation 



Zg _ E::=-ooEr.=-ooexp(-27r^feT^[(m + l/2)^ + n^]) 



(66) 



The KT transition is characterized by Psw — ^/n. For the neighbourhood of 
the KT transition we get 

ZjZp = 0.0864272337... - 0.426489404... (/? - 2/7r) + ... . (67) 

At the KT transition the Psw depends on the scale as 

where in finite size scaling, we identify this scale with the lattice size. It follows for 
lattices with L1 — L2 — L: 

Za/Zp = 0.0864272337... - 0.135755793... / + ... . (69) 

In 1j -\- 



40 



12 Acknowledgements 



I like thank W. Janke for discussions and encouragement. Thanks to E. Vicari for 
reading the final version of the draft. I am grateful to M. Kimball for remarks on 
the status of experiments and for providing me with explicit results for reduced KT 
transition temperatures. This work was supported by the Deutsche Forschungsge- 
meinschaft (DFG) under grant No JA 483/23-1. 



References 

[1] Wilson K G and Kogut J, The renormalization group and the e-expansion, 
1974 Phys. Rep. C 12 75 

[2] Fisher M E, The renormahzation group in the theory of critical behavior, 1974 
Rev. Mod. Phys. 46 597 

[3] Fisher M E, Renormahzation group theory: Its basis and formulation in sta- 
tistical physics, 1998 Rev. Mod. Phys. 70 653 

[4] Pelissetto A and Vicari E, Critical Phenomena and Renormalization- Group 
Theory, 2002 Phys. Rept. 368 549 [cond-mat/0012164| 

[5] Lipa J A, Swanson D R, Nissen J A, Chui TCP and Israelsson U E, Heat 
Capacity and Thermal Relaxation of Bulk Helium very near the X-Point, 1996 
Phys. Rev. Lett. 76 944 

[6] Lipa J A, Nissen J A, Strieker D A, Swanson D R and Chui TCP, Specific 
heat of liquid helium in zero gravity very near the X-point, 2003 Phys. Rev. B 
68 74518 |co nd-mat/0310163j 

[7] Singasaas A and Ahlers G, Universality of static properties near the superfluid 
transition in ^He, 1984 Phys. Rev. B 30 5103 

[8] Goldner L S and Ahlers G, Superfluid fraction of ^He very close to Tx, 1992 
Phys. Rev. B 45 13129 

[9] Barmatz M, Hahn I, Lipa J A, and Duncan R V, Critical phenomena in mi- 
crogravity: Past, present, and future, 2007 Rev. Mod. Phys. 79 1 



41 



Campostrini M, Hasenbusch M, Pelissetto A, and Vicari E, Theoretical esti- 
mates of the critical exponents of the superfluid transition in He4 by lattice 
methods, 2006 Phys. Rev. B 74 144506 |cond-mat/0605083| 

Hasenbusch M, The three-dimensional XY universality class: A high precision 
Monte Carlo estimate of the universal amplitude ratio A^/A^, 2006 J. Stat. 
Mech. P08019 |c ond-mat/0607189] 

Hasenbusch M, A Monte Carlo study of the three-dimensional XY universality 
class: Universal amplitude ratios, [ arXiv: 08 10 . 2 71 6j 

Kosterlitz J M and Thouless D J, Ordering, metastability and phase transitions 
in two-dimensional systems 1973 J. Phys. C 6 1181; Kosterlitz J M, The critical 
properties of the two-dimensional xy model, 1974 J. Phys. C 7 1046 ne- or Two- 
Dimensional 

Jose J V, Kadanoff L P, Kirkpatrick S and Nelson D R, Renormalization, 
vortices, and symmetry-breaking perturbations in the two-dimensional planar 
model, 1977 Phys. Rev. B 16 1217 

Amit D J, Goldschmidt Y Y and Grinstein G, Renormalisation group analysis 
of the phase transition in the 2D Coulomb gas, Sine-Gordon theory and XY 
model, 1980 J. Phys. A 13 585 

M. N. Barber Finite-size Scaling in Phase Transitions and Critical Phenomena, 
Vol. 8, eds. C. Domb and J. L. Lebowitz, (Academic Press, 1983) 

Finite Size Scaling and Numerical Simulation of Statistical Systems, ed. V. 
Privman, (World Scientific, 1990). 

M.E. Fisher, Critical Phenomena, Proceedings of the International School of 
Physics "Enrico Fermi, Varenna, Italy, Course LI, edited by M.S. Green (Aca- 
demic, New York, 1971) 

Capehart T W and Fisher M E, Susceptibility scaling functions for ferromag- 
netic Ising films, 1976 Phys. Rev. B 13 5021 

Gasparini F M, Kimball M O, Mooney K P, and Diaz-Avila M, Finite-size 
scaling of^He at the superfluid transition, 2008 Rev. Mod. Phys. 80 1009 

Dohm V, The superfluid transition in confined 4 He: Renormalization-group 
theory, 1993 PHYSICA SCRIPTA T49A 46. 



42 



[22] Nho K and Manousakis E, Heat-capacity scaling function for confined super- 
fluids, 2003 Phys. Rev. B 68 174503 |cond-ma t/0305500] 

[23] K. Binder, Critical behaviour at surfaces in Phase Transitions and Critical 
Phenomena, Vol. 8, eds. C. Domb and J. L. Lebowitz, (Academic Press, 1983) 
p. 1. 

[24] H. W. Diehl, Field-theoretical Approach to Critical Behaviour at Surfaces 
in Phase Transitions and Critical Phenomena, edited by C. Domb and J.L. 
Lebowitz, Vol. 10 (Academic, London 1986) p. 76. 

[25] Diehl H W, Dietrich S, and Eisenriegler E, Universality, irrelevant surface op- 
erators, and corrections to scaling in systems with free surfaces and defect 
planes, 1983 Phys. Rev. B 27 2937 

[26] Newman K E and Riedel E K, Critical exponents by the scaling-field method: 
The isotropic N-vector model in three dimensions, 1984 Phys. Rev. B 30 6615 

[27] Janke W and Nather K, Monte Carlo simulation of dimensional crossover in 
the XY model, 1993 Phys. Rev. B 48 15807 

[28] Schultka N and Manousakis E, Crossover from Two- to Three-Dimensional 
Behavior in Superfiuids, 1995 Phys. Rev. B 51 11712 |cond-mat/9406014| 

[29] Schultka and E. Manousakis, Scaling of the superfluid density in superfluid 
films, 1996 J. Low. Temp. Phys. 105 3 [cond-mat/96020 85| 

[30] Schultka N and Manousakis E, Boundary effects in superfluid films, 1997 J. 
Low. Temp. Phys. 109 733 [cond-mat/9702216 | 

[31] Zhang C, Nho K, and Landau D P, Finite-size effects on the thermal resistivity 
of^He in the quasi-two-dimensional geometry, 2006 Phys. Rev. B 73 174508 

[32] Hucht A, Thermodynamic Casimir Effect in ^He Films near T^: Monte Carlo 
Results, 2007 Phys. Rev. Lett. 99 185301 jarXiv:0706.3458j 

[33] Vasilyev O, Gambassi A, Maciolek A, and Dietrich S, Monte Carlo sim- 
ulation results for critical Casimir forces, 2007 Europhys. Lett. 80 60009 
[arX iv:0708.2 902j 

[34] Hasenbusch M, The Binder Cumulant at the Kosterlitz-Thouless Transition, 
2008 J. Stat. Mech. P08003 [aFXiv: 0804. 1880] 



43 



Ballesteros H G, Fernandez L A, Martin-Mayor V, and Munoz-Sudupe A, Fi- 
nite size effects on measures of critical exponents in d=3 0(N) models, 1996 



Phys. Lett. B 387 125 [ cond-mat/9 606203| 

Cucchieri A, Engels J, Holtmann S, Mendes T and Schulze T, Universal ampli- 
tude ratios from numerical studies of tie tliree-dimensional 0(2) model, 2002 
J. Phys. A 35 6517 [cond-mat/02020lf| 

Deng Y, Blote H W J, and Nightingale M P, Surface and bulk tran- 
sitions in three-dimensional 0(n) models, 2005 Phys. Rev. E 72 016128 
[cond-mat/0504173i 

Bittner E and Janke W, Nature of Phase Transitions in a Generalized Complex 
\psi\^ Model, 2005 Phys. Rev. B 71 024512 |cond-mat/0501468| 

Hasenbusch M and Torok T, High precision Monte Carlo study of the 3D 
XY-universality class, 1999 J. Phys. A 32 6361 |cond-mat/9904408| 

Campostrini M, Hasenbusch M, Pehssetto A, Rossi P, and Vicari E, Critical 
behavior of the three-dimensional XY universality class, 2001 Phys. Rev. B 63 
214503 | cond-mat/0 010360i 

Li Y-H and Teitel S, Finite-size scaling study of the three-dimensional classical 
XY model, 1989 Phys. Rev. B 40 9122 

Fisher M E, Barber M N, and Jasnow D, Helicity Modulus, Superfluidity, and 
Scaling in Isotropic Systems, 1973 Phys. Rev. A 8 1111 

Hasenbusch M, Direct Monte Carlo Measurement of the Surface Tension in 
Ising Models, 1993 J. Phys. I (France) 3 753 | hep-lat/92090l 6| 

Wolff U, Collective Monte Carlo Updating for Spin Systems, 1989 Phys. Rev. 
Lett. 62 361 

Hasenbusch M, Pinn K and Vinti S, Critical Exponents of the 3D Ising Uni- 
versality Class From Finite Size Scaling With Standard and Improved Actions, 
1999 Phys. Rev. B 59 11471 |h ep-lat/9806012j 

Saito M, An Application of Finite Field: Design and Implementation of 128-bit 
Instruction- Based Fast Pseudorandom Number Generator, PhD thesis, Dept. 



44 



of Math., Graduate School of Science, Hiroshima University, Advisor: M. Mat- 
sumoto; The numerical program and a detailed description can be found at 



'|http: / /www.math.sci.hiroshima-u.ac.jp/~m-mat/MT / SFMT / index.html ' 



Weber H and Minnhagen P, Monte Carlo determination of the critical temper- 
ature for the two-dimensional XY model, 1988 Phys. Rev. B 37 5986 

Hasenbusch M, The two-dimensional XY model at the transition temperature: 
A high precision Monte Carlo study, 2005 J. Phys. A: Math. Gen. 38 5869 



|cond-m at/0502556l 

Hasenbusch M and Pinn K, Computing The Roughening Transition of Ising 
and Solid-on-Solid Models by BCSOS Model Matching, 1997 J. Phys. A 30 



(1997) 63 (cond-mat/9605019] 



Caselle M and Hasenbusch M, Universal Amplitude Ratios in the 3-D Ising 
Model, 1997 J. Phys. A 30 4963 [higUit/9701007| 

Sabisky E S and Anderson C H, Onset for SuperRuid flow in He^ Films on a 
Variety of Substrates, 1973 Phys. Rev. Lett. 30 1122 

van de Laar R W A, van der Hoek A, van Beelen H, Two-dimensional behaviour 
of very thin ^He films on glass, 1995 Physica B 216 24. 

Yu Y Y, Finotello D and Gasparini F M, Finite-size scaling and the convec- 
tive conductance and specific heat of planar helium films near the superfluid 
transition, 1989 Phys. Rev. B 39 6519 

Kimball M O, private communication 

Gasparini F M and Metha S Temperature oszillations from a temperature 
superleak, 1998 J. Low Temp. Phys. 110 293 

[56] Kimball M O and Gasparini F M, Superfluid fraction of He-3-He-4 mixtures 
confined at 0.0483 fi m between silicon wafers, 2001 Phys. Rev. Lett. 86 1558 



45 



